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Abstract 


Unified constitutive material models were 
developed for structural analyses of aircraft gas 
turbine engine components with particular applica- 
tion to isotrrpic materials used for high-pressure 
stage turbine blades and vanes. Forms or combina- 
tions of models independently proposed by Bodner 
and Walker were considered in this study. These 
theories combine time-dependent and time- 
inoependent aspects of inelasticity into a continu- 
ous spectrum of behavior. This is in sharp con- 
trast to previous classical approaches that parti- 
tion inelastic strain into uncoupled plastic and 
creep components. Predicted stress-strain 
responses from these models were evaluated against 
monotonic and cyclic test results for uniaxial 
specimens of two cast nickel-base alloys, B1900+Hf 
and Reng 80. Tension-torsion test results for 
Ren§ 80 and Hastelloy X alloys were used to evalu- 
ate multiaxial stress-strain cycle predictions. 

The unified models, as well as appropriate algo- 
rithms for integrating the constitutive equations, 
were implemented in finite-element computer codes. 

Introduction 


This paper discusses the development of uni- 
fied constitutive material models and of compatible 
finite-element computer codes for the structural 
analysis of gas turbine engine hot section compo- 
nents. This effort constitutes a different 
approach for nonlinear structural analysis which 
has heretofore been based on classical plasticity 
theory with uncoupled creep equations. Previous 
analytical studies of hot section components such 
as turbine blades' and combustor liners^ have 
demonstrated that existing nonlinear finite-element 
computer codes baseo on classical methods do not 
accurately predict the cyclic response of the 
structure. 


The unified constitutive theories considered 
in this study encompass time-dependent and time- 
independent aspects of inelasticity including 
plasticity, creep and stress relaxation. These 
theories avoid the noninteractive summation of 
inelastic strain into plastic and creep components 
or specifying yield surfaces to partition stress 
space into elastic and elastic-pla tic regions. 

In eliminating these overly simpliiied assumptions 
of classical theory, unified models can more 
realistically represent the behavior of materials 


under cyclic loading conditions and high tempera- 
ture environments. 

Results are reported for the first 18 months 
of a 4 yr NASA-sponsored program to develop inelas- 
tic structural analysis methods based on unified 
constitutive theories. This work is being per- 
formed under two NASA contracts; one with the 
General Electric Company and the other a joint 
effort of Southwest Research Institute and Pratt & 
Whitney. Model development is directed toward the 
prediction of the structural response of initially 
isotropic, cast nickel-base alloys for the temper- 
atures and strain ranges characteristic of air- 
cooled turbine blades and vanes in advanced gas 
turbine engines. 

After extensive review and evaluation of uni- 
fied theories which have been suggested in the 
literature, those proposed by Bodner 3 and Walker^ 
were selected for further study and development. 
Both of these theories were significantly revised 
by their authors during the course of and as a 
result of this program. Forms or combinations of 
these two models, as well as numerical methods for 
integrating the constitutive equations, have been 
implemented in two finite-element computer codes; 
by G.E. in a modular form in a new code specifi- 
cally developed for a unified theory and by PtW in 
a stand-alone module for use with general purpose 
codes such as MARC. 

Both programs followed essentially the same 
experimental work plan. This consisted of uni- 
axial and multiaxial tests to establish data 
bases, develop efficient procedures for determin- 
ing the material constants and evaluate the candi- 
date models. The alloys tested were Ren£ 80 in 
the G.E. program and B1900+Hf in the SwRI/PfcW pro- 
gram. More than 30 uniaxial, isothermal, strain- 
controlled tests were performed under various 
strain rates, temperatures, alternating-to-mean 
strain ratios and hold times. The tests were 
continued until cyclic saturation was achieved. 

In adoition, three uniaxial thermomechanical 
cyclic tests were conducted to study the models 
under nonlsothermal conditions. Most of the 
multiaxial tests involved tension-torsion testing 
although G.E. also performed some cyclic tests on 
notched bar specimens to simulate tension-tension 
stress fielos. Comparisons of model multiaxial 
predictions are shown for Ren6 80 and for previ- 
ously conducted tension-torsion test results for 


Hastelloy X alloy. As a final demonstration of 
the unified analytical methods, cyclic structural 
analyses will be performed for an aircooled tur- 
bine blade or vane airfoil under simulated engine 
operating conditions. 

Constitutive Models 

The Initial effort In this program was to 
survey and evaluate unified constitutive theories 
proposed in the literature. Only those theories 
were considered that Incorporated time- Independent 
(plasticity) and time-dependent (creep) aspects of 
inelastic behavior In a single Inelastic strain 
'•ate term. Also excluded from consideration were 
theories which required a yield surface. 

The criteria used In screening candidate 
constitutive theories were (1) ability to 
accurately represent the cyclic inelastic behavior 
of an isotropic, high-temperature superalloy, 

(2) applicability over the temperature and strain 
ranges of Interest for aircooled turbine blades 
and vanes and typical mission cycles of advanced 
airbreathing engines, (3) capability of extension 
to multlaxial stress states and anisotropic materi- 
als, (4) reasonableness of material property test- 
ing requirements and ease of determining model 
constants, and (5) integrability of constitutive 
equations for implementation In finite-element 
computer codes. 

Most unified models can be described by a set 
of constitutive equations that have the basic form 

(S - ft) 

M) 

k - h, (K)ij - r, (2) 

U = 1^(0 " ^ tft )^! 1 ” r 2 ( 3 ) 

Equation (1) is a flow law relating the 
inelajtic strain rate and the stresses where S 
and c. are deviatoric stress and Inelastic strain 
rate tensors. The tensor internal variable, o , 
defines the kinematic or directional hardening 
(the Bauschinger effect) and is frequently refer- 
red to as a back stress. K is a scalar internal 
variable, commonly called a drag stress, which 
defines the isotropic hardening. Temperature 
effects on the inelastic strain rate are generally 
taken into account by considering some of the 
material constants of the constitutive model to be 
temperature dependent. Equations (2) and (3) are 
evolutionary equations describing the growth laws 
for the internal variables. Most evolutionary 
equations for the back stress-drag stress type of 
model include both hardening and recovery terms 
where h] and h 2 are functions describing 
the hardening, d is a dynamic recovery function 
and r] and r 2 are static thermal recovery 
functions. Almost all back stress-drag stress 
models use the inelastic strain rate, £,, as the 
hardening criterion. About a dozen unified theo- 
ries which have been proposed in the literature 
were considered in the survey. There were sub- 
stantial differences among the models in the func- 
tional relationships used in the constitutive 
equations. 


The major exception to the basic form of 
Eqs. (1) to (3) is the Bodner model. The flow law 
for this model is of the form 



where 0 and n are material constants with 0 
representing the limiting strain rate in shear 
ana = 1/2 

A major Difference between Eqs. (1) and (4) is 
that the back stress models assume the direction of 
the inelastic strain rate vector to be coincident with 
the direction of (£ - q) whereas the Booner model 
assumes it to be coincident with direction of £. 

Isotropic and directional hardening are dis- 
tinguished by a partitioning of the internal vari- 
able Z into components Zj and Zo rather 
than by a back stress. 

Z - Zj ♦ Z D (5) 

The evolutionary equations for the internal 
variable are of the form 

*i = h i (z i - Z A - r i l6 > 

* D = w fi P - d( vS, - r 2 {1) 

The constant Z] in Eq. (6) is the saturation 
value of Zj. Cyclic hardening or softening is 
controlled by the Zj component and depends on 
whether the initial value of Zj is less than or 
greater than Zj. The procedure for incorporating 
directional hardening in the theory was developed by 
Bodner during this program. Another difference 
between the evolutionary equations for the Bodner 
model and the back stress-drag stress type of mgdel 
is that the former uses the plastic work rate, W p , 
as the measure of hardening rather than the inelas- 
tic strain rate. There are essentially nine mate- 
rial constants to be determined for the Bodner 
model, of which only three were founo to be tem- 
perature-dependent . 

Because of its unique form and its simplicity 
compared to the more common back stress models, 
the Bodner model was selected for more detailed 
study under both contractual efforts. Of the back 
stress-drag stress type of theory, both contrac- 
tors selected the Walker model as a basis for fur- 
ther study. The reason for this selection was 
that the Walker model was further along in develop- 
ment for finite-element analysis than other com- 
parable back stress-drag stress models. The major 
difficulty in applying the Walker model is to 
determine the 14 temperature-dependent material 
constants. 

In the course of this study, 8odner and Walker 
(who were used as consultants in the SwRI/PWA pro- 
gram) incorporated additional terms in their models 
to account for nonproportional loading. As the 
measure of nonproportionality, Bodner used the 
angle between the stress and stress rate directions 
while Walker used the angle between the strain and 
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strain rate directions. Walker also revised his 
moael to use an exponential law In place of the 
previous power law for the functional form of the 
kinetic term, f(ii), in Eq. (1). 

A composite of the two models was developed 
In the G.E. program with the purpose of retaining 
the basic simplicity of the Bodner model and using 
the directional hardening form of the Walker 
model. This Involved incorporating a back stress 
In the flow law of the Bodner model to give 



where k 2 * 1/2(S - o)(S - a) and A Is a constant 
related to the afctivfctiTin energy. The back stress 
growth law for the composite model, as shown below, 
is essentially Eq. (3) with an additional function, 
g, to Improve the modeling of the cycle In the 
elastic range. 

r * Wtl - “Mill - r Z 

£ *[i - gluty* ♦ g(fl)$] (io) 

The directional form of Z given by Eq. (7) was 
eliminated. Isotropic hardening Is described by 
the Isotropic form of the Z variable given by 
the evolutionary Eq. (6). There are six material 
constants to be experimentally determined for the 
composite model, neglecting thermal recovery terms. 

evaluation of Constitutive Models 

Eva’uation of the candidate unified theories 
was based on extensive uniaxial and multlaxlal 
testing. The parameters varied in these tests 
were temperature, strain rate, strain range, mean 
strain, dwell time and, for the multlaxlal tests, 
stress field. Details of the test specimens, con- 
ditions, and procedures are described in the first 
annual reports 5 *” of the contractors. One of the 
significant differences between the two contract 
efforts was 1r the cyclic behavior of the base 
materials that were selected. The G.E. alloy, 

Ren4 60, exhibited some cyclic softening at ele- 
vated temperatures while the SwRI/P&W alloy, 
B1900+Hf, underwent cyclic hardening. The discus- 
sion below is directed mainly toward the Bodner 
model for the SwRI program and the composite model 
for the G.E. program. Difficulties were encoun- 
tered in both programs in determining the Walker 
model constants. 

Bodner Model 

All of the constants for the B 1 900+Hf alloy 
were obtained from uniaxial, monotonic, tensile 
stress-strain tests. These tests were conducted 
over a range of temperatures from room temperature 
up to 1093 # C and over a range of strain rates 
from 0.0005 to 0.1 min"'. A ccmoarison of the 
calculated and experimental stress-strain curves 
is shown in Fig. 1 at temperatures of 538, 871, 
and 982 # C for a strain rate of 0.005 min* 1 . To 
determine the constants associated with strain 
hardening, the tensile data were correlated in 
terms of the derivative of the stress with respect 
to the plastic work plotted against stress, as 


shown In Fig. 2. The curves for each strain rate 
are seen to be essentially parallel. Each of 
these curves was approximated by two linear 
segments. From the slopes and Intercepts of the 
segments, the Isotropic (Zj) and directional ( Zq) 
variables could be readily determined. The slopes 
will change with temperature at the higher temper- 
atures due to thermal recovery effects (Eqs. (6) 
and (7)). For a cyclic softening material, this 
correlation would result In entirely linear curves 
and some cyclic tests would be required to deter- 
mine the directional component (Zq). 

Cyclic response was successfully predicted 
for the B 1 900+Hf alloy using the material con- 
stants obtained from the monotonic tensile stress- 
strain tests. As shown in Fig. 3, the cyclic 
predictions during initial cycling and at 
saturation show excellent qualitative agreement 
with the experimental cycles. 

An essential requirement for a unified theory 
is the capacity to characterize material creep 
behavior. This is demonstrated for the Bodner 
model in Fig. 4 where predicted and experimental 
results are compared for uniaxial, constant-load 
creep tests of 81 900+Hf at various temperatures. 

As with the cyclic prediction, presented in 
Fig. 3, the material constants were derived from 
monotonic tensile stress-strain tests at different 
strain rates. 

The ability of the Bodner model to account 
for cyclic hardening under nonproportional loading 
is shown In Fig. 5. Comparison of the experi- 
mental (Fig. 5(a)) with the predicted (Fig. 5(b)) 
cyclic response of Haste! loy X tension-torsion 
specimens subjected to 90° out-of -phase loading 
indicate good correlation. The Bodner model 
always assumes coincidence between the deviatorlc 
stress and plastic strain rate directions as shown 
in Fig. 6 while the tension-torsion test cycles 
had phase angles from 7° to 45°. However, this 
discrepancy in phase angle between the model 
assumptions and multiaxial experiments did not 
appear to adversely affect the accuracy of the 
cyclic predictions. 

Walker Model 

The material constants of the Walker model 
were less aeterministic ana more sensitive to the 
strain rate than the constants for the Booner 
model. The moael could predict the cyclic 
response for kcn£ BO adequately provided a dif- 
ferent set of constants was determined for several 
strain rates. This limitation in representing a 
range of strain rates for Ken£ 80 cyclic behavior 
applied to similar back stress-orag stress models 
that were investigated, but not to the composite 
model . 

Cyclic stress-strain predictions for 
Hastelloy X under 90° ojt-of -phase tension-torsion 
loading (Fig. 6(c)) show good agreement with the 
experimental results. Although the predicted 
phase angles in Fig. 6 are 5° to 15° smaller than 
the experimental observations, this is more real- 
istic than the assumption of coaxiality between 
the deviatorlc stress and plastic strain rate vec- 
tors for the Bodner model. However, the theoreti- 
cal results shown in Fig. 5(c) and Fig. 6 were 
obtained by fitting the Walker model to the exper- 
imental response in Fig. 5(a), whereas the Bodner 
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model predictions were based only on uniaxial 
monotonic test data In determining the material 
constants. 

Composite Model 

To better represent the cyclic softening 
behavior of Ren£ 80, G.E. developed a composite 
model which retained the functional form of the 
Bodner model but Incorporated a Walker-type back 
stress In the constitutive equations. The growth 
law for this back stress was based on the Inelas- 
tic strain rate as in the Walker theory while 
retaining the Bodner Isotropic haraening variable 
based on the work hardening rate. All the mate- 
rial constants, except one, were determined from 
monotonic tensile stress-strain tests. The con- 
stant controlling the cyclic softening behavior 
was determined from saturated cyclic hysteresis 
loops. Figure 7 compares the calculated and ex- 
perimental ken£ BO *♦: ss-strain curves at 382 ®C 
for strain rates of u.002, 0.02, 0.06, and 
0.2 min -1 . The agreement is seen to be very 
good, especially considering the data scatter 
indicated by the two sets of tests at a strain 
rate of 0.2 min -1 . 

Composite theory predictions of the high- 
temperature cyclic behavior of Reng 80 are com- 
pared with experimental hysteresis loops in Fig. 8 
for fully reversed cycles at 982 °C. The stress- 
strain loops in Fig. 8(b) are for a relatively 
slow strain rate of 0.002 min -1 and in Fig. 8(a) 
for a high strain rate of 0.2 min -1 . The model 
predictions were able to capture the cyclic 
response reasonably well for both strain rates. 

The composite model was also evaluated with 
respect to its ability to predict the creep char- 
acteristics of the material. In Fig. 9 inelastic 
strain calculations are compared to test data for 
constant load creep tests at 982 °C. As can be 
seen, the steady-state creep rate was reasonably 
predicted, but the primary creep rate was too 
high. The creep analysis capability of the model 
was also assessed for a strain controlled cycle 
involving stress relaxation, as shown in Fig. 10. 
In this test case, a 120 sec hold time was applied 
at the minimum strain limit of the cycle. Compar- 
ison of the analytical and experimental hysteresis 
loops show that the calculateo minimum stress part 
of the cycle was too low, and as a consequence the 
calculated stress relaxation was too rapid. How- 
ever, since the analytical and experimental loops 
were in agreement at the end of the stress relaxa- 
tion, it appears that the model worked reasonably 
well in predicting the stress relaxation. 

The composite model has been used to predict 
the response of Rend 80 at 9B2 °C under propor- 
tional and nonproportional multiaxial loading con- 
ditions. Figure 11(a) shows the stress behavior 
under a 90 degree out-of -phase tension-torsion 
cycling. It is seen that the model has the abil- 
ity to analyze complex cyclic hardening and soft- 
ening behavior under nonproportional loading 
conditions. Figure 11(b) shows the predicted 
phase angle between the inelastic strain rate and 
the deviatoric stress vectors. Comparison of 
Figs. 6 ano 11 indicates that the phase angles are 
similar although the tests were very different. 
However, the loading conditions in the two multi- 
axial tests were similar and it has been found 
that the phase angle is a strong function of the 


loading condition. Further refinement of the 
composite model is expecteo as more multiaxial 
test data Become available. 

Implementation of Models in Finite-Element Codes 

Unified constitutive theories involve mathe- 
matically stiff differential equations where small 
changes in the independent variables can cause 
large changes in the dependent variables. Success- 
ful implementation of unified theories in finite- 
element codes requires the use of efficient methods 
for integrating the constitutive equations. The 
suitability of any integration method is dependent 
or. the iteration and incrementation algorithms 
used in the computer code ano on the nature of the 
problems to be analyzed. In both contractual 
efforts an initial strain technique was employed 
for derating the equilibrium equations. Since 
the constitutive flow law and the evolutionary 
equations are uncoupled, a further iteration is 
requireo within the overall equilibrium iteration 
to solve for the inelastic strain rate and inter- 
nal variables. The iteration process causes the 
solution to tend to drift from the equilibrium 
state. To correct for this drift, residual load 
corrections are applied to the incremental loads 
by computing the difference between the external 
and internal forces summed over all the elements 
at the end of the previous increment. Automatic 
time incrementation has also been incorporated in 
the incrementation procedures to allow for time 
step control under large cyclic stress and inelas- 
tic strain excursions. The control is achieved by 
limiting the maximum stress and inelastic strain 
change permitted for any time increment. 

General Electric developed both two- and 
three-dimensional nonlinear finite-element codes. 
These codes are modular and incorporated different 
inelastic constitutive models. A dynamic time 
incrementation procedure is used to minimize 
cost. The constitutive equations are integrated 
using a secona order Aaams-Moulton predictor- 
corrector technique. The number of iterations 
required has been reduced ty improving the stabil- 
ity of the initial strain procedure. Figure 12 
shows the analysis of an Inconel 718 benchmark 
notch specimen' using the Booner model with the 
two-oimensional code. It can be seen that the 
cooe nut only performed well, but providea excel- 
lent agreement wch the measured notch root 
strains. 

In the SwRI/P&W program, the constitutive 
equations for both the Bodner and revised Walker 
models, as well as the computational algorithms, 
were incorporated in the MARC user subroutine, 
Hypela. With some modifications this subroutine 
can become a stand-alone module for use with other 
general purpose nonlinear codes. A number of 
methods for integrating the Bodner and Walker 
constitutive equations are still being studied. 

An explicit Euler technique with a self-adaptive 
solution method has been found to be suitable for 
both models. 

Concluding Remarks 

Tne continuing oevelopment and revisions of 
unified constitutive material theories, as well as 
differences in the cyclic behavior of the materi- 
als stuoieo makes difficult any Definitive judge- 
ment at this point in the program as to which 
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constitutes the best approach. However, the 
results up to the present time give encouragement 
that this program will meet Its objective to sig- 
nificantly advance the structural analysis capa- 
bilities for engine hot section components. These 
Improvements In the structural analysis tools 
snoula provide better life prediction capability 
and greater hot section durability for the next 
generation of gas turbine engines. 
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Figure 4. - Results from constant load creep 
tests for 81900 +Hf. 
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Figure 5. - A comparison of Bodner and Walker theories with experimental results on the stress-strain 
response of Hastelloy-X tested under 90 P out-of-phase tens ion /tors ion loading at ror m temperature. 




Figure 6. - A comparison of Bodner and Walker theories 
with experimental results on the phase angle between 
the deviatoric stress and the plastic strain rate vectors 
observed in Hastelloy-X tested under 90 P cut-of-phase 
tension/torsion loading at room temperature. 
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Figure 7. - A comparison of the calculated (composite theory) and 
the experimental stress-strain curves of Rene 80 at 982 °C. 
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Figure 8. - A comparison of the calculated 
(composite theory) and the experimental 
hysteresis loops for Rene 80 at 982 °C. 
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Figure 9. - Results from constant load creep tests of 
Rene 80 at 982 °C. 



Figure 10. - A comparison of the calculated 
(composite theory) and the experimental 
hysteresis loops with 120 sec. hold time 
at minimum strain for Rene 80 at 962 °C. 


PHASE ANGLE, deg SHEAR STRESS 12. MPa 


ORIGMNAl 

OF POOR QUALITY 



-300 -200 -100 0 100 200 300 

STRESS 11, MPa 


(a) Cyclic softening behavior. 



(b) Phase angle between the deviatonic stress and inelastic 
strain rate vectors. 

Figure 11. - Composite model predictions of Rene 80 at 982 °C 
under 90P out-of-phase tens ion /tors ion loading. 
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Figure 12. - Comparison of predicted and measured notch root strains using the Bodner model with 2 D computer 
code developed at General Electric. 
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